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We derive the equations for the gravity assist manoeuvre in the general 2D case without the constraints of 
circular planetary orbits or widely different masses as assumed by Broucke 1 , and obtain the slingshot conditions 
and maximum energy gain for arbitrary mass ratios of two colliding rigid bodies. Using the geometric view de- 
veloped in an earlier paper by the authors 2 the possible trajectories are computed for both attractive or repulsive 
interactions yielding a further insight on the slingshot mechanics and its parametrization 3 . The general slingshot 
manoeuvre for arbitrary masses is explained as a particular case of the possible outcomes of attractive or repul- 
sive binary collisions, and the correlation between asymptotic information and orbital parameters is obtained in 
general. 



I. INTRODUCTION 

The slingshot or gravity assist manoeuvre—"— is often consid- 
ered as part of a restricted three-body problem and its use 
has been associated in the literature mostly with spaceflight 
strategies^"— with some applications in astrophysics for the 
study of mass ejection from binary clusters 11 and the proposal 
of new General Relativity testsi^. In reality, the design of 
spacecraft trajectories between two planets is a many-body 
problem except for the slingshot part, which is in most de- 
signs well approximated by an elastic binary collision. This 
work focuses on the slingshot manoeuvre as a particular case 
of a general binary elastic collision between massive objects 
subjected to central interaction forces. In a previous work 2 the 
geometric determination of binary collisions was introduced, 
and the possible outcomes were in some cases surprising 3 . We 
have obtained a parametrization of all possible outcomes of 
a binary elastic collision in an arbitrary frame, and from the 
mass ratios and initial velocities as asymptotic initial condi- 
tions we obtain a picture not only of the final asymptotic ve- 
locities in terms of a single parameter 6 in the 2D case, but also 
the detailed description of the two-body motion that fits these 
asymptotic data in the case of the gravitational or coulombian 
interaction. We are therefore in condition to determine which 
precise orbital parameters must be chosen to obtain a desired 
effect on a flyby of a satellite about a planet or star, be it a 
gravity assisted boost or capture. The conditions for a gravity- 
assisted manoeuvre of a satellite are often loosely associated 
to a flyby in front or behind the planet. 1 Both in the case of 
attractive as well as repulsive collision forces, a harder look 
must be performed to really grasp what the critical ingredient 
is. In particular we show how the geometry and the timing of 
arrival at the point of closest approach (periapsis), to wit the 
relative position of the bodies with respect to the normal to the 
Center of Mass velocity V cm at that point, determines the out- 
come of the collision and refine the phenomenological rule-of- 
thumb that a flyby in front of the planet results in a slowing of 
the satellite whilst a flyby behind the planet's trajectory would 
result in a boost. 13 In fact, the asymptotic description of a col- 
lision is somewhat elusive in this respect. By scaling out the 



interaction in all space and time dimensions there is naturally 
a loss of information of precisely where and when are the two 
bodies for given initial velocities, so a determination must be 
made as to what corresponds in a real problem to these times 
and positions. It could be claimed that these initial velocities 
should be those of the bodies when their distance is equal to 
the sum of the radii of their respective spheres of influence 14 . 
But since this too is a fuzzy concept this is not much of an 
improvement. In fact, that information is only present when 
enough conditions are specified to determine the collision out- 
come uniquely. 15 Thus, in the 2-dimensional case, the circum- 
ference of possibilities for the velocity outcomes of one of 
the bodies encodes the missing information about where the 
bodies initially are simultaneously when they have the given 
velocities. This can also be translated into an impact param- 
eter in the non-inertial body-frame for one of the masses (or 
reduced-mass frame) but that begs the question of viewing the 
collision in the laboratory frame. In this paper the assumption 
will therefore be made that at t = the Center of Mass (CM) 
will be at the origin of the laboratory frame (LF). 

Diagrams like the one in Fig. Q] introduced in a previous 
work£, are used to correlate the asymptotic information, pro- 
vided by initial velocities far away from the periapsis, with 
the possible eccentricities, focal distances and other orbital 
parameters for open Keplerian orbits in case of gravitational 
attractive or Coulombian repulsive scattering. These diagrams 
depict the relation between incoming laboratory frame asymp- 
totic velocities v D and u , of masses ra v and m u respectively, 
and their final asymptotic velocities Vj and u p through a com- 
putation involving the scattering angle 6 of the m u mass, mea- 
sured in its initial asymptotic M-body frame from the direc- 
tion of the incoming relative velocity v Q - u of the m v mass. 
In that M-body frame the circumferences of possible veloc- 
ity outcomes are easily drawn and their image in the labora- 
tory frame can easily be deduced, thus yielding information 
about the possible directions and magnitudes of asymptotic 
final velocities Vj and iij for both masses. In this way the 
orbits can be viewed in the laboratory frame and a study can 
be made, for instance, for the optimal incidence angle on a 
planetary fly-by that delivers the maximum allowed velocity 
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boost in a chosen direction. The energy gain is obtained both 
in the case of idealized point particle collisions and extended 
object collisions, where the periapsis distance is constrained 
by a minimum value below which the collision is no longer 
elastic. The relation of this asymptotic information with the 
actual trajectory can also be displayed using the fact that the 
reduced-mass-frame trajectory axis and asymptotic directions 
are already included in these diagrams, and thus also all the or- 
bital parameters can be deduced or introduced here. The final 
velocities can be determined from these diagrams by specify- 
ing the scattering angle 8 of u l - u relative to the inbound 
direction v D - u in the initial frame of one of the bodies (the 
M-body S u frame in Fig. [TJ. This angle specifies the asymp- 
totic direction of the scattering of the w-body in the S u frame 
where it was initially at rest, the reference direction being de- 
fined by the CM velocity V cm - u = ^;(v - u Q ) in that 
frame. Denoting by ti fl the direction defined by d; — u , and 
w fl its orthogonal direction in the plane of the collision, we 
deduce from the momentum and kinetic energy conservation 
laws that 2 
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together with the orthogonality equations 
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In the interesting limiting case where to v » to u and v D « 
V cm » Vj we can anticipate that the \ x circumference reduces 
to a point at V cm , meaning that the to v body motion is prac- 
tically unaltered by the collision. Then d ^ D = 2 (v - u ) 
and, for given magnitudes v , u D , this is maximized for head- 
on collisions which provides the greatest magnitude variabil- 
ity for the outbound u 1 . This is the case of planetary flyby 
by satellites for attractive orbits, and it is immediately appar- 
ent that scattering in the direction of the planet's velocity (if 
possible) has the potential for more dramatic acceleration or 
deceleration of the satellite. This is however limited by the 
fact that planets have large radii (eventually including an at- 
mosphere) and the satellite cannot get closer at periapsis than 
that radius. Thus some values for 8 may be excluded as un- 
practical. 



II. DETERMINATION OF ORBITAL PARAMETERS 



The binary elastic collision diagrams such as that shown in 
Fig. [Uprovide several simultaneous views of the event, to wit 
the laboratory, Center of Mass and body frames. The pos- 
sible velocity outcomes for both masses are parametrized by 
the angle 8 that d[ - u D makes with the reference direction 
v Q - u , which is also the direction of the CM velocity in the 
S u initial body frame. In this frame, where m u was initially 
at rest, 8 is the scattering angle of the mass to u after the colli- 
sion with an incoming mass to v . The range of 8 e [— |, +| ] 
encompasses all possible results in a binary elastic collision 
with given initial velocities. Thus 8 works in the S u frame 
as the extra parameter needed to determine the outcome of 
the collision uniquely, a role that is usually attributed to the 
'impact parameter' b (the distance in reduced-mass-frame be- 
tween the origin and the asymptotic line drawn from the in- 
coming body, with direction this body's asymptotic relative 
velocity). This impact parameter can only be precisely de- 
fined in the reduced-mass frame, which we will assume here 
to be the instantaneous body-frame S v of the mass to v (since 
we are ultimately interested in exploring all directions for the 
fly-by about a planet of mass to v > to u , the satellite will be 
henceforth represented by the M-body of mass m u ). 

In the S Y body-frame (BF) with axes parallel to those 
of the laboratory frame (LF), the relative motion of the 
two masses will appear as that of a single reduced mass 
/j. = to v to u /(to v + to u ), at the relative position f " of to u 



(6) 



which appears to be moving under interaction forces pointing 
to a fixed total mass M = m v + m u at the origin, where the 
other mass (to v ) is at rest, and this can actually be computed 
for sufficiently well behaved forces. 

Once the (BF) motion v'^(t) is obtained, and assuming that the 
frame directions are parallel to those of the laboratory frame, 
we may return to the (LF) description by noting that, in the ab- 
sence of external forces, the CM motion is uniform and there- 
fore 

*v<0 = R ™(o - = R cm (0) + V cm f - ^ *;<*), 



f„(o = R cm (o + ^f;w = R cm (0) + v cm ? + ^ ?;(?). 



(7) 

This is an approximation in the real case of planetary flyby be- 
cause of the gravitational influence of the sun, but for the du- 
ration of the encounter, assumed to start and end at the bound- 
aries of the planetary sphere of influence, the effect of the third 
body is assumed to be negligible. 16 
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FIG. 1: Collision diagram for different masses m v > m lx with initial velocities \ , u , and its relation with the collision diagram as seen from 
the w-body (mass m u ) initial rest frame S u . In this frame the circumference with diameter d = 2(V cm - u Q ) is the locus of all possible outcomes 
Uj - u parametrized by the scattering angle 6. A choice of 6 determines u l - u and an orthogonal vector ii restricted to a circumference 
of diameter D = 2(v D - u Q ). The outbound velocity \ t - u is geometrically defined as ii + ii, — v and lies on a circumference centered on 
V cra - u with radius \\ a - V cm |. The angles 6 and ip represent the u- and v-scattering angles relative to the incoming velocity of the CM in the 
K-body initial rest frame and iff is the scattering angle in the CM frame. 



A. View from the v-body frame S v (reduced mass 

SYSTEM) 

If the (BF) reference directions were rotated and scaled with 
respect to the (LF) directions in such a way that r = A ■ f ", 
with AeOj the rotation matrix, then we should write (0 in 
the form 

f v (o = R cm (0) + v crn f-^A.f;(o, 

(8) 

_i _^ III 

?„(o = R cra (0) + v cm f + ^A.?;(o. 

Viewed from the reduced-mass frame S Y attached to the v- 
body of mass m v with reference directions rotated through L 
by A" 1 ((/>), the asymptotic incoming velocity u^' of the w-body 



of mass m u obeys (see Fig. [2]i 

A-u" = U =u - v , (9) 

o o o o ' v 7 

which is anti-parallel to the diameter vectors d and D of the 
reference circumferences. Likewise, the asymptotic outgoing 
velocity u" for the w-body in this frame verifies 

A ■ u" = Uj = u 1 - Vj . (10) 

On the other hand, in the reduced-mass frame S v the mean 
force acting on the M-body during the collision is proportional 
to that body's total linear momentum variation Ap^' in said 
frame and 

a • Ap; = ft a • (a," - a;) = ^ (a, - u + v - Vl ) . (ii) 

Since conservation of linear momentum in the laboratory 
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frame implies 

Av = v, -v D = -H u -u, = --pAu, (12) 

Ml ' 1-1-7 x ' 



(fTTb is always the sum of collinear terms parallel to iij - u G , 
and in fact 



• AP^' = j« (Au - Av) = m u (fij 



Ap„ 



(13) 



In the reduced-mass frame S v , for Newtonian or Coulombian 
type interactions, the asymptotic motion is known to be conic 
with one focus at the origin 17 . In particular for strictly positive 
total initial energy £>" > in the (BF), a gravitational colli- 
sion must be an hyperbola concave to the focus at the origin, 
while for repulsive Coulombian forces the hyperbola branch 
is convex. The direction of Ap (( in the (LF) then also deter- 
mines the hyperbolic axis e" for the M-body trajectory in S v , 
and its angle i) w with the asymptotes -u^' and u" is also 
related to the parameter angle 8 (see Figs. [2] and [3). Once 
the angle (f> between the axis e A " oc Ap^' of the hyperbola and 
a chosen laboratory frame e r direction is known, the rotation 
matrix A((h) such that r = A • r" is determined as 



I cos(0) 
sin(0) 




- sin(0) 
cos(0) 
1 



(14) 



assuming that the constant angular momentum is L = L_e_. 

According to the diagrams in Fig. [2] for a given incoming 
direction u " in the S v frame, making an angle 



4> = tan 



(15) 



with the e v direction in the LF frame, the directions (f> = O + ( 
of the possible axis t" oc Ap " are in the range 
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(16) 



The sign of <p - <p D = 6 is an indication of whether the incom- 
ing motion is from above or below the axis 



(17) 



The angle <f> = O corresponds to the situation where e" is 
aligned with u^', i.e. a head-on collision with 6 = 0. 



In the S v frame, the angle &_ m from the e " axis to the asymp- 
totic direction -u" = A"'(0) • U Q is d_ m = tan" 1 (-u^') = 6 
and it determines the boundaries for the actual orbit. For 
a repulsive interaction, the polar angle <? r = & will change 
in the interval & T e while f° r an attractive in- 

teraction & a = n — d and its domain of variation is 
# a e [ji - 0_ m , 2n - as indicated in Fig. [3] 



B. Initial conditions and determination of orbit 

PARAMETERS 



Knowing the asymptotic angle d = n - of the hyperbolic 
trajectory f^'(f) = f "(f) of the reduced mass fi in the v-body 
frame S v determines the eccentricity e of this orbit through 
the relation 
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(18) 



The actual branch of the hyperbola that corresponds to the 
motion is the one concave towards the focus "Fj at the origin 
O if the interaction is attractive, otherwise it is the one concave 
to the other focus !F 2 , situated along the axis e" at a distance 
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(19) 



from the hyperbolic center. Here s is the focal distance (to 
the directrix) and a the semi-major axis (i.e. half the distance 
between axis intercepts with the hyperbola) see Fig. [3] 

Using polar coordinates \r" , 8\ in the S v body frame, an hy- 
perbola with axis e A " aligned with A _1 (0) • e t and a focus at 
the origin O = f\ will be defined for all eccentricities e > 1 
and focal distances s > by the parametric equation 



£ S 

^ ^ ~ 1 + e cos(y ± c?) 



(20) 



where y is a constant dependent of initial conditions. Us- 
ing henceforth k = G m v m u for a gravitational interaction, or 

k = | | for a Coulombian interaction, the time equation can 
also be expressed as a parametric function of # in the form of 
Kepler's equation 



2 tanh 
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e sin(#) 
1 ± e cos(#) 



(21) 



where the sign choice distinguishes attractive (+) or repulsive (-) orbits. 
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FIG. 2: Asymptotic directions and axis (blue) as viewed in the v-body frame. This is a M-body slingshot collision where m v «s 2.6m u . Notice 
the interpretation of u, - u o as proportional to Ap„, the total momentum change of the M-body in the reduced-mass S r frame, and also as the 
symmetry axis of the trajectory in the non-inertial v-body frame. The angle then is identified with the asymptotic angles d ±oo of the trajectory 
in that frame. 



From the eccentricity e and the initial energy in the S v frame, 
£" = \ fi liig'l 2 , it is possible to obtain the focal distance s and 
the angular momentum £. " 



k (e 2 - 1) _ k sin(6») tan(0) 



"|2 



"|2 



(22) 



£' 



/^(e 2 -l) k 2 tan 2 (0) 



;"i2 



i"|2 



and consequently (since £" = n b |u "|) the impact parameter 
b and the displacement a e of each focus from the origin. 



b = 
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III. THE SLINGSHOT MANOEUVRE 

In Fig. [2] we have an example of a collision where the lighter 
body gains kinetic energy as seen from the laboratory frame. 
This is a near-maximum slingshot collision for M-body with 
mass m u jk 0.3m v for given initial asymptotic conditions. This 
is assuming unrestricted periapsis conditions, i.e. point parti- 
cle collision. In real collisions not all 8 angles are accessible 
in the vicinity of jift for the outgoing M-body because that 
would imply a periapsis distance smaller than allowed by the 
dimensions of the bodies for an elastic collision. 

The slingshot manoeuvres are particular cases of the possible 
outcomes for either an attracting or repelling collision. It is 
possible to realize them in all collisions with arbitrary mass 
ratios, even though it only provides significant boosts in cases 
where a small inertial mass collides with a much larger one 
moving much faster. The nature of the interaction is irrele- 
vant, as long as it is central and conservative. Evidently we 
are not considering for the moment variations of the manoeu- 
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FIG. 3: Orbital detail of the attracting and repelling open Keplerian 
or Coulombian orbits viewed in the v-body frame <S V . The e" axis is 
chosen to coincide with the reduced mass total change-of-momentum 
Ap^' in the collision. The asymptotic angle 6 is the same as the one 
identified in figs[TJand[2] 



vre such as aero-gravity assisted slingshots 18 where at perigee 
the forces involved are neither conservative nor central. The 
particular value of 8 that corresponds to the theoretical max- 
imum slingshot case is when 6 equals half the angle if/ D be- 
tween v Q - u D and V cm , which means that fi[ and Vj would 
both come out collinear to the CM velocity V cm . This can be 
shown as follows. From one of the orthogonality conditions 
in equations (0 written as 



llfi] -u ll 2 = d-(d! ~u ) 



with 



d = 2fv -u)= 2Mv (y -u ), 



(24) 



(25) 



one gets after expanding on both sides, 

"i 2 " 2uj • u D + u\ = 2V cm • (d, - u ) - 2u, • u D + 2u 2 Q . (26) 

therefore 

Ul 2 = u 2 +2V cm \\u l -uJIcosOA), (27) 

where if/ is the angle from iij - u to V cm . But from the first 
equation (01 



2m„ 
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llv -u o llcos(0), 



(28) 



where is the angle from v D - u D to iij - u Q . Then, if i// 
denotes the angle from v o - u D to V cm , the relation 

<A = -Ao - 6 (29) 

holds for if/ Q g [-ji, n] and 8 g [-§ , |]. Then the magnitude 
of U[ depends solely on 6 and ip Q through the relations d27T i to 



u\ = u 2 + V cm \\y a - u ll cos(<A - 8) cos(8) . (30) 

m„ + m„ 



From ( f30b we conclude that the range of 8 that provides for a 
M-body velocity boost or slingshot (corresponding to u \ — > 
0) is bounded by boost-break angle 
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which is defined through the boost condition 
cos(i/r - 8) cos(#) > 0. 



(31) 



(32) 



This means that m a -boosts will happen for 8 g [ - | , ^ bb [ 
if if/ < (and for 8 g ] 8 hb , | ] if > 0).ln each case, the 
remainder of the ^-domain will correspond to a breakage of 
the w-body and a boost to the v-body velocity. From the defi- 
nition of 8 hb we conclude that for i(j o = + all collision results 
are m u -velocity boosts (for instance when m u < m v these are 
"head-on" collision with opposing velocities |u | < |v | and 
arbitrary impact parameter), while for tf/ D = ±n all scenarios 
correspond to a m u -velocity break (when m u < m y this is a col- 
lision where m u "catches-up" m w with parallel velocities and 
arbitrary impact parameter). Not surprisingly, when v - u D 
is perpendicular to V cm we have equal 8 domains for getting a 
boost or a breakage. 

The condition d32b above also means that cos(i/r) > 0, that 
is, in any slingshot situation we find that the angle tp from 
fi[ - u Q to V cm must verify \tf/\ < |. Since Ap„ = m u (u l - 
u D ), besides representing the total impulse acting on the u- 
body during the maneuver, is also the direction of the force 
acting on it at the periapsis, the slingshot condition can now 
be phrased as follows: if at the point of closest encounter the 
force acting on the M-body has a positive component in the 
direction of the CM velocity, then there will be a boost in 
the final M-body velocity. Otherwise we will obtain a M-body 
velocity break. Notice that this formulation is valid for both 
gravity assisted and coulombian slingshots, i.e. attractive as 
well as repulsive interactions. 

The angle if/ can be obtained from initial conditions in terms 
of the v to u D angle f3 Q (see Fig. 0, in which case we can 
express 8 bb by 



tan 



(1 + V)Xo sin^o) 



(33) 



where rj = — and^- = ^. Notice that when rj = 1 and^- Q = 1 
then tf/ = | for all fi , and 8 bb = 0, which means that in every 
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collision we have equal 8 domains for boosting or breaking. 
For the more common planet-satellite case 77 » and Eqn. 
( l33T > reduces to 



0, 



bb 



tan 



cscQgJ 



(34) 



The following figures depict three typical situations in gravity 
assisted manoeuvres. These were derived from the equations 
deduced so far for the binary collisions and using diagrams 
such as those in Figure (fj). A live Java applet that models 
these collisions with a variety of mass ratios and zooming 
scales can be found in [— ]. The first two figures correspond 
to slingshot boosts of the lighter body m u , while the last one 
corresponds to a breaking manoeuvre of m u . Figure|4]-(a) rep- 
resents a slingshot boost for bodies with similar velocities. 
Figure[4]-(b) is a "catch-up" collision where the heavier, faster 
body boosts the lighter body velocity as it passes by it. Fig- 
ure 3]-(c) is a breaking collision for similar velocity bodies. 
Besides the asymptotic velocities they show the actual trajec- 
tories near the periapsis. The color coded points represent 
same-time position in both orbits. Adjacent point intervals do 
not correspond to equal time intervals but rather equal angular 
displacements in the reduced-mass frame orbital representa- 
tion. The zero index point corresponds to the periapsis. As 
can be seen from these examples the usual rule-of-thumb for 
boost or break according to wether the lighter mass passes be- 
hind or in front of the heavier one can be very tricky to apply 
when the mass ratio is not too big. Our proposed formulation 
above in this section is however rigorous and unambiguous. 



But this then means that the outgoing asymptotic direction 
dj — Vj coincides with V cm . Denoting these slingshot extremes 
by uf and Vj 1 then there is a scalar A such that 



S W Sl mi„ 

ui ' - vi = AV r 



(39) 



Since in general \\u l — v x II = llu - vjl, the previous equation 
yields 



A = m °~^ 



I IV,. 



(40) 



From the definition of V cm the following holds also in general 



u, - V, 



-(fii-V cm ), (41) 



so in conjunction with (f40b 
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and finally the maximum slingshot velocity is 
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m llu — V " ^ 

\ _|_ V Q O 



(42) 



(43) 



cm / 



m v + m u ||V c 

Corresponding to this maximum, the velocity 
vi = ui - AV cm will be the minimum of possible 



A. Slingshot maxima 



For ideal point masses the unconstrained extremes are found 
as usual through 



dl (u\) = u\ d\\j + d g u\ d8 = (d g u\ - u\) d8 = 0, 



meaning 

sin((/r) cos(6) - cos(i/0 sin(0) = sin(<^ - 8) = . 



(35) 



(36) 



This holds if if/ - 8 = nn for integer n, but physically only 
n = and n = 1 are of interest. There are thus two extremes. 
The second variation of u\ defines their type depending on 
whether & 2 (u\) ^ 0. But 

d 2 (uj) = d sin(</r - 8) = -2 cos(t// -8)d8, (37) 

so if/ - 8 = corresponds to a maximum while iff - 8 = n 
corresponds to a minimum. 

Thus the maximum slingshot boost for given incoming initial 
conditions will happen when fij - u cc Ap 2 makes an equal 
angle if/ = 8 with both V cm and v - u . According to ( 1291 this 
is 



1 



1 



cos 



(38) 



sl m 
Vl 



(44) 



cm / 



On the other hand when if/ = 8 + n then V cm points in the 
opposite direction of iij - Vj, so instead 



Vl - Ul 



and then 



Vl 



^sl nlin 
Ul 



1 + m u ""o-V ' 

m v +m u \\\ j 



1 - 



* cm ' 



m v + m u \\\ 



(45) 



(46) 



(47) 



IV. THE CONSTRAINED SLINGSHOT 



The preceding calculations are however only valid when the 
periapsis distance of the hyperbolic orbit can be taken as small 
as necessary, a situation that does not apply in most physical 
applications. In fact, usually the intervening bodies have di- 
mensions which prevents an elastic collision to occur if the 
periapsis distance r p is smaller than a limiting value r p < r min . 
This can be for instance the radius of a planetary mass or 
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(c) 

FIG. 4: Orbital detail of trajectories for: (a) an attractive interaction near a maximum slingshot of a M-body (black) of mass m u colliding 
elastically with a mass m v ss 3m u . (b) attractive slingshot of a M-body (black) of mass m u in a "catch-up" collision with a mass m v « 2m u . 
(c) breaking maneuver for M-body of mass m u « |m v . Color coded dots represent simultaneous (same color) positions in the trajectory, in 
equally-spaced ^-intervals. The closest approach corresponds to label 0. 



star, plus an arbitrary offset to prevent friction from eventual 
atmospheres or stellar corona. Thus a different maximiza- 
tion must be performed which embodies this constraint when 
calculating the maximum kinetic energy boost per unit mass 

l 2 9 

Ak H = j(mi - Ug) that can be achieved for orbits that do not 
exceed the limit of maximum approach given by a specific 



We will first do this calculation in the limiting case of V cm « 
v which is typical for gravity assist flyby of a small satellite 
about a planetary mass, i.e. when m u <k m y . 
Recall that for an hyperbolic orbit the periapsis distance is the 
r-value parametrized by & = 0, that is 



Tp l+e' 



The semi-major axis a and the focal distance s 



a = 



L" 
Kfie ' 



(48) 



(49) 



can be used to express the kinetic energy £" = i U% = k/2 a 
and r p = (e - I) a from which follow the identities 



a 



a = ^T„ (50) 

where we recall k = Gm w m u . For an attractive orbit 
cos((?_ oo ) = -l/e really means cos (6) = 1/e (since d_ m = n-6) 
that is 



cos(#) 



1 



1 



1 + 



1 + 



(51) 



Now, since V cm » v we obtain the following identities relat- 
ing initial velocities (see Fig. [5]): 



^o sin (^ ) = M o sin (A,)> 
cos(iAJ = v Q -u Q cos{fi ) . 



(52) 



where /3 a denotes the incidence angle of u relative to the 
direction v Q and \f/ is the angle between the directions v and 
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U D = u D - v D . We can use d52l ) to eliminate U a = llv -u ll and 
tp from Eqn. (l30b and express the change in kinetic energy 
per unit mass of the satellite in the laboratory (heliocentric) 
frame as 

A£ B = cos(0) 2 [l - i (cosOSJ - sin08 o )tan(e))j . (53) 



Using Eqn. dBTT l to eliminate 8 from Eqn. d53l ) we can express 
the gain in kinetic energy per unit mass of the satellite in terms 
of the approach angle /3 and the velocity ratio x Q = U J V - 



2»< 



1 -Xo cos 08 o ) - J (l + $ ixl - 2 cosOSJ^ + l)) - 1 
. * i 

( 1 + ^fc 2 -2cos^ o)/ro + l)) 2 



(54) 



Notice that in ( f54b under the assumed conditions [i/m u « 1 
and «/r p jU » G is the potential energy per unit mass at 
the periapsis distance of the planet. 

Now we will derive the general relation for the energy gain per 
unit mass in situations where the momenta of the bodies are 
of similar order in magnitude, in which case V cm is no longer 
collinear with v„. In this case the angular relations presented 
in Fig. [5]indicate that the convenient decomposition for \p is 



4> = <A - 6 = a + OAo - a ) ■ 



(55) 



where (t/r o - a) is the angle between V cm and v Q . Furthermore, 
using the definition of V cm , the projection V cm • v D can be 
expressed as a function of /3 D and the ratio x a = M </ v o 

^cm V o cos («Ao - a) = H v l | — + —Xo cos^o) 1 - (56) 



Now we should write (see Figure [5]) 

( U sin(a) = u Q sia(fi ), 

I U Q cos(aO = v D - u Q cosQSJ , 

{ U^v^tf -2cos08 o ) Ao + l) 



(58) 



Substitution of Eqns. d55l l to d58b and (T5Tb into formula (1301) 
gives for equal masses m v = m u = m 



Ak„ = 



( 1 + ^ta 2 -2* cos(/y + i)) 



(59) 



and in the general case 



and also we can get from \ a x V cm that 



V cm v sin(<A -a) = -^v 2 oXo sin(^ ) . (57) 
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Ak. = 



2»X 



Ju 1 " Xo + 2Xo sin08 o )^(i + ^-(^^^cos08 o ) + i))^i 
(l + ^fc 2 -2^cos08 o ) + l)) 2 





( 1+ ^rU 2 - 2 ^ cos ^ )+1 ) 


2 \ 

-i 




'l + ^(^o-2A-oCos08 o ) + l)) 2 



(60) 



Notice that in these last formulas we have assumed that the 
periapsis distance is fixed at r p . What this distance is depends 
entirely on the particular nature of the collision type. For a 
slingshot problem of a small satellite about a planet r p might 
be a small multiple of the planet radius, but for a binary star 
collision it should probably be greater than the Roche limit to 
avoid the inelastic effects of the gravitational tides between 
massive extended bodies. 

Formulas (|54| |. (|59l ) and ( l60t are generalizations to arbitrary 
masses and orbital conditions of the known energy change ex- 
pressions for gravitational slingshot such as those found in 
Broucke's paper 1 . 



V. APPLICATION TO GRAVITY-ASSIST MANOEUVRES 



In the particular case of spacecraft manoeuvres assisted by the 
gravitational field of a planetary object we can assume that 
m v » m u in the formulas above, therefore to a high accuracy 
^- « 1 and - K » gi— . When approaching a large planetary 

mass M = m v from an angle j3 a (measured from v D to u , see 
Fig. [3), a satellite of mass m = m u cannot engage the collision 



with an arbitrary periapsis distance r p and has to maintain a 
minimum distance larger than a multiple of the planet's ra- 
dius R. Thus the extreme limiting periapsis distance will be 
assumed as the planetary radius, r p > R. From the expression 
for the inverse eccentricity - = cos(#) we obtain a limitation 
on the available orbits by specifying that possible collision 
outcomes must respect the relation 

OTS(e)< i + ^iiLvir (61) 

CM o T o" 

or, in terms of the ratio x a = 7 1 of satellite to planetary speeds 

1 



cos(#) < 



1 + ^(1-2^)^+^) 



(62) 



This relation limits the maximum velocity outcome at the 
slingshot through Eqn. [53] otherwise it could theoretically 
be 



sl m 
Ml 



v a (l + A /l-2cos(/3 ) A - 0+Ao 2 j (63) 

if there were no limitations on the periapsis distance allowed. 
In terms of relative kinetic energy gain per unit mass we can 



1 1 > 




FIG. 6: Boost scenario for a collision with a mas- 
sive tiny object for which ^ « 1. 



express this for the theoretical case of an encounter with a 
massive object with very small radius as 



2Ak„ 



= 1 



2\2 



Xo 



(64) 



for which we can expect the general boost pattern as in Figure 
[6] For solar system planets however the situation is not as 
favorable, and we get instead the scenarios shown in Figure [7] 
below. Notice that in both cases the conditions for a breaking 
manoeuvre are much more limited than those for the boost. 
For other planets the scenarios are similar to these two cases. 
These manoeuvres are purely gravitational and do not include 
the effects of powered boosts at periapsis to modify orbital 
parameters and thus achieve a different orbit, or aerogravity- 
assisted manoeuvres to alter the bending angle and yield larger 
boostsJ2 
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Earth (Ak„xlO — ) 
Kg 



Jupiter (Ak„xlO — ) 
Kg 





FIG. 7: Energy gain per unit mass for Earth and Jupiter gravity-assist manoeuvres for limiting-case trajectories with periapsis equal to the 
planetary radius. 



VI. CONCLUSIONS 



We have shown that the slingshot effect is no stranger than any 
other elastic collision, just a particular case of the general set 
of possible outcomes where the timing of arrival at the point 
of closest approach is more favorable to a boost in the veloc- 
ity of the lighter mass. In the case of point particles, equation 
(|43T > indicates that the optimal policy is a quasi-head-on for 
attractive interactions, and a head-on collision for repelling 
interactions. However, for extended bodies there is a mini- 
mum distance r min of approach beyond which there is a severe 
departure from elasticity (e.g. crash and burn), therefore the 
approach angle 8 for u - v Q must provide for sufficient ec- 
centricity e such that (a - l)e > r min for attractive ones, or 
(a + 1) e > r min for repulsive interactions. 

In the gravitational case, the Hohmann transfer orbits 20 may 
be the more energetically efficient but they are not necessarily 
the optimum policies for approaching a slingshot configura- 
tion. Still in this case there is often the question of whether 
there could be a solar slingshot manoeuvre. The answer is 
obviously affirmative. In fact, just like any other elastic bi- 
nary collision, a slingshot around the Sun is a possibility when 
viewed from a frame where the Sun itself is moving. Also 
obvious is nonetheless the fact that, just as in the case of 
the Jupiter gravity-assist fly-by where there is no gain in the 
satellite velocity when viewed from the planetary frame, any 
slingshot manoeuvre involving the sun would appear to return 
a disappointingly unaltered final speed in the solar reference 
frame, but this need not be the view from another frame. How- 
ever it is true that for most velocities u' we are able to throw 



spacecraft with towards the Sun, as seen from that external 
frame where the Sun itself is moving with velocity v « V cm 
the spacecraft velocity would be 



u = u' + V 

o o cm 

and thus the boost term in (T43b would be 
m„ IliL-vJI llu 



m„ + m„ 



(65) 



(66) 



So unless llu^ll is already much greater than the speed of the 
solar system itself, this will in general be a small factor and 
the resulting speed of the spacecraft is about the same as that 
of the Sun itself. There is also the question of producing a 
positive energy orbit for a spacecraft that is already within the 
sphere of influence of the Sun without which the incoming 
orbit would not be truly hyperbolic. 

Note that angular momentum has no role in these calculations, 
other than being a globally conserved constant that defines the 
plane where the collision takes place and defines the relation 
between the impact parameter b and the parameter 6. So it 
is surprising to see that many so-called didactic explanations 
of the effect still mention 'stealing angular momentum' to ac- 
count for the increase in velocity of a spacecraft in a gravity 
assisted fly-by. It should also be stressed that, although in the 
specific case of planetary fly-by the calculation of actual or- 
bits involve a complex A^-body problem, there is no aspect of 
the slingshot effect above that involve more than two-body in- 
teractions. In particular no three-body effects are needed to 
understand the slingshot of spacecraft in the vicinity of plan- 
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etary masses 21 , although the calculations for intermediate tra- 
jectories and timings for launch and arrival at the planet vicin- 
ity should include this for added accuracy. Tisserand's cri- 
teria and three-body graphical methods have been proposed 
for these high-precision calculations 22 , but the fact remains 
that the slingshot effect involves basically the mechanics of 
an elastic binary collision. 

Having said this, mention should also be made to the limi- 
tations on the presented model for the calculations of an ac- 
tual planetary fly-by. The provisos made in [ 2 ] concerning the 
asymptotic nature of the collision process must be pondered 
with actual data: the times involved in the approach and ex- 
traction of a spacecraft from the fly-by should be compared 



to the duration of the fly-by before they can be considered as 
asymptotically infinite. For the duration of the fly-by, the ex- 
ternal forces must be negligible as compared to the collisional 
interaction to introduce only minor perturbations in the result- 
ing trajectories, since it takes too long to disregard the fact 
that both the planet and the spacecraft are orbiting the Sun. 
Thirdly, the entry in the planet's 'sphere of influence' (which 
is often considered to be the Hill's sphere) depends on the ap- 
proach that is made, and that should also be factored in the 
calculations. Still, as long as these factors can be shown to in- 
troduce small perturbations to the simplified binary collision 
model, this can be used to successfully explain the physics of 
the effect. 



* Electronic address: amaro@fisica.ist.utl.pt 
T Electronic address: lemos@fisica.ist.utl.pt 

1 R. A. Broucke. The celestial mechanics of gravity assist. In Amer- 
ican Institute of Aeronautics and Astronautics, editor, AIAA/AAS 
Astrodynamics Conference, Minneapolis, MN, Aug. 15-17, 1988, 
number (A88-50352 21-13) in Technical Papers, pages 69-78, 
August 1988. 

2 A. Rica da Silva and J. P. S. Lemos. Geometric parametrization 
of binary elastic collisions. American Journal of Physics, 74(7): 
584-590, 2006. 

3 Amaro Rica da Silva. 2D Elastic Collision 
Parametrization (An Interactive Java Applet). 
http://centra.ist.utl.pt/ amaro/Collisions/Collisions.html , 2005. 

4 John J. Dykla, Robert Cacioppo, and Asim Gangopadhyaya. 
Gravitational Slingshot. American Journal of Physics, 72(5):619- 
621, May 2004. 

5 Kenneth J. Epstein. Shortcut to the Slingshot Effect. American 
Journal of Physics, 73(4):362, April 2005. 

6 J. M. Longuski and S. N. Williams. Automated design of gravity- 
assist trajectories to mars and the outer planets. Celestial Mechan- 
ics and Dynamical Astronomy, 52:207-220, 1991. 

7 Robert C. Johnson. The slingshot effect. 
http://maths.dur.ac.uk/ dmaOrcj/Psling/sling.pdf January 2003. 

8 G. D. Racca. New challenges to trajectory design by the use of 
electric propulsion and other new means of wandering in the solar 
system. Celestial Mechanics and Dynamical Astronomy, 85: 1-24, 
January 2003. 

9 V. V. Malyshev, V. E. Usachov, and Y. D. Tychinskii. Solar Probe 
Mission with Multiple Gravity-Assist Maneuvers Realized with 
Conversion Launchers. Cosmic Research, 41(5):431 - 442, 2003. 

10 C. A. Ocampo. Transfers to earth centered orbits via lunar gravity 
assist. Acta Astronautica, 52(2003): 173 - 179, 2003. 

11 W. C. Saslaw, M. J. Valtonen, and S. J. Aarseth. The gravitational 
slingshot and the structure of extragalactic radio sources. Astro- 



physical Journal, 190:253-270, 1974. 

12 J. M. Longuski, E. Fischbach, and D. J. Scheeres. Deflection of 
Spacecraft Trajectories as a New Test of General Relativity. Phys- 
ical Review Letters, 86(14), April 2001. 

13 A. V. Labunsky, O. V. Papkov, and K. G. Sukhanov. Multiple 
gravity assist interplanetary trajectories, volume 2 of Earth Space 
Institute book series on public and private sector interest in space. 
TF-CRC, 1988. ISBN 90-5699090-X. 

14 V. Barger and M. Olsson. Classical Mechanics: A modern per- 
spective. McGraw-Hill, 1995. 

15 H. Asada. An exact solution to determination of an open orbit. 
Celestial Mechanics and Dynamical Astronomy, 97:151 - 164, 
March 2007. 

16 E. Barrabes, G. Gomez, and J. Rodriguez- 
Canabal. Advanced topics in astrodynamics. 
http://www.ieec.fcr.es/astro04/notes/gravity.pdf, Barcelona, 
July 2004. 

17 A.E.Roy. Orbital Motion. Institute of Physics Publ., 2005. 

18 R. Armellin, M. Lavagna, and A. Ercoli-Finzi. Aero-gravity assist 
maneuvers: controlled dynamics modeling and optimization. Ce- 
lestial Mechanics and Dynamical Astronomy , 95:391 -405, May 
2006. 

19 W. R. Johnson and J. M. Longuski. Design of Aerogravity-Assist 
Trajectories. Journal of Spacecraft and Rockets, 39(l):23-30, 
2002. 

20 William E. Wiesel. Spaceflight Dynamics. McGraw-Hill, 1997. 

21 James A. Van Allen. Gravitational assist in celestial mechanics-a 
tutorial. American Journal of Physics, 71(5):448-451, May 2003. 

22 J. K. Miller and Connie J. Weeks. Application of Tisserand's cri- 
terion to the design of gravity assist trajectories. In AIAA/AAS 
Astrodynamics Specialist Conference and Exhibit, AIAA 2002- 
4717, August 2002. 



